Entanglement entropy of two disjoint blocks in critical Ising models 
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We study the scaling of the Renyi and entanglement entropy of two disjoint blocks of critical Ising 
models, as function of their sizes and separations. We present analytic results based on conformal 
field theory that are quantitatively checked in numerical simulations of both the quantum spin chain 
and the classical two dimensional Ising model. Theoretical results match the ones obtained from 
numerical simulations only after taking properly into account the corrections induced by the finite 
length of the blocks to their leading scaling behavior. 
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Conformal field theory (CFT) is one of the most pow- 
erful and elegant tools to study quantum one-dimensional 
(ID) systems and classical two-dimensional (2D) ones. It 
provides a complete description of the low-energy (large- 
distance) physics of critical systems that can be classi- 
fied only on the base of their symmetries [1]. One spec- 
tacular recent success was the application of this frame- 
work to 2D turbolence [2]. The predictions of CFT have 
been tested in experiments for carbon nanotubes [3] , spin 
chains [4], and cold atomic gases [5], just to cite a few of 
the most recent ones. 

CFT has been traditionally applied to the computation 
of large distance correlations of local observables. Only 
recently it has been realized that CFT is also the ideal 
tool to describe the global properties of a large subset 
of microscopical constituents (e.g. spins) and in partic- 
ular their entanglement. This has generated an enor- 
mous interest in the study of the entanglement proper- 
ties of many-body systems [6] that is connecting several 
branches of physics such as quantum information, con- 
densed matter, black hole physics. The quantum infor- 
mation insight about the origin of the achievements of the 
density matrix renormalization group (DMRG) in ID, 
and its failure in higher dimensions [7], can be cited as 
an example of the outstanding results generated by this 
cross-over between different branches of physics. The en- 
tanglement between two complementary regions A and 
B of a quantum system described by the state \i[>) can 
be measured through the entanglement entropy. This is 
defined as the von-Neumann entropy of the reduced den- 
sity matrix p A = Tr B \yj) obtained by tracing over the 
degrees of freedom in the region B. When vb is the ground 
state of an infinite ID critical system and A is a block of 
length £, CFT predicts the universal scaling [7-9] 

S A = C -\ogl + c' 1 , (1) 

where c is the central charge and c[ a non universal con- 
stant. This formula is the most effective way to calculate 
the main signature of the CFT (the central charge), and 
it can be used to identify the universality class of new 
models, as for example done in the Fibonacci chain [10]. 



The reason of this simple scaling in CFT is easily un- 
derstood [9]. In fact, through a replica trick, Sa can be 
interpreted as — d n Trp A \ n =i. For integer n, Trp^ is the 
partition function on an n-sheeted Riemann surface with 
two branch points at the border of the interval A that can 
be mapped to the plane by a conformal transformation. 
By studying the transformation of the stress-energy ten- 
sor under this conformal mapping, one has that Trp^ is 
the two-point correlation function of some twist- operators 
that have scaling dimension A n = c/24(n — 1/n), i.e. 
Trp'X = c n £- c / e ( n - 1 / n '> . By analytically continuing this 
to complex n and by taking the derivative in 1, we get 
Eq. (1). This reasoning also applies to the case of TV 
intervals: Trp A is the partition function of a n-sheeted 
Riemann surface with 2N branch points, i.e. a 2iV-point 
function of the same twist-operators. A generally incor- 
rect result was obtained by uniformizing this surface [9]. 
This is not allowed because of the non-zero genus of the 
Riemann surface. This result was checked in several free- 
fcrmionic theories [11] and only recently, the error has 
been pointed out [12-14]. In the case of many inter- 
vals, Trp A turns out to be a function of the full oper- 
ator content of the theory and not only of the central 
charge. For a free compactified boson or Luttinger liq- 
uid (LL) Trp A has been calculated for n = 2 [13] and 
for general integer n [14]. However, the functional de- 
pendence on n is so complicated that the analytic con- 
tinuation has not yet been achieved. These predictions 
have been checked against the exact diagonalization of 
the XXZ chain [13, 14]. Unfortunately, the numerical 
results are limited to relatively small system sizes and 
only few general properties (like the dependence on the 
LL parameter) have been checked: large oscillating cor- 
rections to the scaling (as for one block [15]), have made 
impossible a quantitative comparison for the scaling func- 
tions related to Trp^. Concepts and calculation schemes 
used to get these results (such as higher genus Riemann 
surfaces, twist fields, orbifold theories) are mathemati- 
cal tools that have been mainly used in string theory 
and that only now find their place in condensed matter 
physics. 



The entanglement of many intervals thus depends on 
the details of the CFT and should be calculated case by 
case [16]. The simplest and most studied CFT is the 
critical Ising model that in the continuum is a free Ma- 
jorana fermion and has central charge c = 1/2. The 
corresponding ID quantum spin chain is the Ising model 
in transverse field described by the Hamiltonian 



(2) 



are Pauli matrices acting on the spin at site j 



where a 

and we use periodic boundary conditions. The model has 
a quantum critical point at h = 1. The correspondence 
with a free fermion could erroneously lead to the conclu- 
sion that Sa for the Ising chain is the incorrect result of 
Ref. [9], valid for free fermion theories [11]. This is not 
the case when the block A involves more than one inter- 
val since the unitary transformation that maps the spin 
degrees of freedom to the fermionic ones is not anymore 
contained inside A, as it is easily checked by direct cal- 
culation [17]. Sa for two intervals has been calculated in 
the Ising chain [18], but for the fermion degrees of free- 
dom and it agrees with Ref. [9]. The breaking of the 
equivalence of fermions and spins makes any lattice ex- 
act computation hard, and a representation of pa for two 
blocks is not yet known. For this reason, we analyze the 
problem with numerical methods. We use a tree tensor 
network (TTN) algorithm [19] for the quantum ID Ising 
model [20] and MonteCarlo simulations of the classical 
2D one as in Ref. [12]. Using the mapping to the torus 
partition function for n = 2, we provide the CFT pre- 
diction for Trp^. The generalization of this result to all 
integer n requires a more detailed analysis (as for the LL 
[14], but more difficult because of the complexity of the 
target space [21, 22]) that we are currently studying and 
will be reported elsewhere [23]. 

We consider the case of two disjoint intervals A = 
[u\ , u 2 ] U [us , U4] . By global conformal invariance Tvp\ 
can always be written as [13, 14] 



F n (x) , (3) 



where u; 



Uj and x 



M21U43 / (U31U42) is the four- 



point ratio. F n (x) is the universal scaling function that 
depends on the theory, and c„ the non-universal factor 
of the single block. The normalization is F n (0) = 1. 
The incorrect result of Ref. [9] is F n (x) — 1 identically. 
For a chain of finite length L, one replaces Uij by the 
chord distance L/n sm(nuij/L). F n (x) is symmetric for 
x -> 1 -x [13]. 

The TTN (as the better known DMRG) gives the full 
spectrum of the reduced density matrix. From this Sa 
and the moments of pa can be extracted and analyzed. 
The scaling functions F n (x) (for the entropy Fvn{x) = 
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FIG. 1: TTN scaling function F^ix) vs the conformal ratio x 
for different block sizes I, The upper points are the extrap- 
olation to £ — > 00 from Eq. (5). Data for L 7^ 512 are not 
shown in the legend. The full line is the CFT prediction (6). 



—F[(x)) are obtained as ratios (difference) of Trp^ (Sa) 
with the prefactor in Eq. (3). We consider two blocks of 
length I at distance r. The four-point ratio x is obtained 
by substituting in its definition the chord distance: 



sin tt£/L 



sin tt( 



(4) 



In the x variable, we would expect that data with dif- 
ferent £, r and L would collapse onto a single curve thus 
revealing the scaling functions F n (x). 

We start our analysis from the data for the function 
F 2 (x) reported in Fig. 1 for I between 2 and 128 and 
L from 64 to 512. The finite I results do not display 
the symmetry x — > 1 — x and the data present large 
corrections to their leading scaling behavior. To extract 
the asymptotic behavior we perform a finite-size analysis. 
For any x, general RG arguments give the scaling 



F^(x,£) = F. 



CFT 
2 



(x)+t- s 'f 2 (x) + 



(5) 



where S c is an unknown exponent, f 2 (x) is the scaling 
function of the first sub-leading correction, and the dots 
indicate further ones. The data are well described by 
5 C = 1/2. The evidence of this scaling for different x is 
shown in Fig. 2. It is easy to extrapolate to I — > 00 
(the points where the straight lines cross the vertical 
axis) and the results are reported in Fig. 1. The ex- 
trapolation restores the symmetry x — > 1 — x. It is pos- 
sible to calculate this quantity from CFT. In fact, the 
2-sheeted Riemann surface has the topology of the torus, 
on which it can be mapped by a conformal transforma- 
tion. The torus partition function for the Ising model 
is 2Zl ms = {p u=2 \9 v {r)/r,{r)\) 2 [1], where r, (r) is the 
Dedekin function, #„(t) are the Jacobi elliptic functions 
and t is the modular parameter. In our case, r is given 
by the solution of x = [9 2 (t) /0 3 (t)} 4 [21]. For this value 
of r, major simplifications occur (as for rj = 1/2 in the 
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FIG. 2: Corrections to the scaling for F 2 at (a;) at fixed x in 
Eq. 5. Inset: universality of fa (x). The dashed line is oc x 1 ^ 4 . 

LL [13]) and the final result can be written in terms of 
only algebraic functions: 



F2(x) = 



1 



1/2 



+ x 1/A + ((1 - x)xY /i + (1 - x) 



1/4 



a/4 



1/2 



(6) 



This curve is reported in Fig. 1 and agrees with incredible 
precision with the extrapolated data. For i « 1 we have 



F 2 (x) = 1 



,1/4 



/2 



In the inset of Fig. 2, we 



report the universal correction to the scaling function 
f 2 (x) obtained as [i^ at (x, t) - F§ FT (x^l 1 / 2 for different 
I that collapse (without any adjustable parameter) on a 
single curve. In the inset we show f 2 (x) ~ re 1 / 4 . 

To check the universality, we study the classical critical 
2D Ising model, using the algorithm of Caraglio-Gliozzi 
to obtain the two-point function of twist- fields [12]. We 
use an asymmetrical geometry with the temporal direc- 
tion Lt equal to 10 times the spatial one L (between 24 
and 324). The results for F 2 {x) are reported in Fig. 3 
showing the same qualitative features as Fig. 1. The 
extrapolations to £ — > oo present large error bars, but in 
agreement with CFT. This also implies that a rescaling 
of all (large enough) length scales should give the same 
numbers in the two models (as in 2D [24]). The rescaling 
factor a can be calculated from the single block entan- 
glement obtaining L 2 d = o,Lid, with a ~ 0.71. In the 
inset of Fig. 3, the MonteCarlo data for the L = 8 clas- 
sical systems are compared with the L = 6(~ 8 x 0.71) 
quantum chain showing a good agreement. 

In Fig. 4 we report the TTN scaling function for 
Fyx{x). Unfortunately the CFT value is unknown be- 
cause we are not yet able to make the analytic continu- 
ation (as for the LL). One important feature is evident 



FIG. 3: MonteCarlo determination of the scaling function 
F 2 (x). The full line is the CFT prediction (6). Inset: Com- 
parison between MonteCarlo data for the 2D classical Ising 
model and the exact diagonalization of the quantum chain. 



from the plot: the corrections to the scaling are negligible 
and all data collapse in a single symmetric scaling curve. 
In the inset of the figure we report the data in log-log 
scale to emphasize the power-law behavior for small x. 
In the LL, F n (x) for small x displays a power-law with 
an n-independent exponent [14]. This reasoning general- 
izes to the Ising model [23] and from the result for F 2 (x) 
we read that the exponent is 1/4, as confirmed by the 
plot. We also found that the prefactor is it. Moreover, 
for various n, we computed the function F n (x) for the n- 
th moment of pa also showing large finite I corrections. 
The analysis of these data will be reported elsewhere [23]. 
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FIG. 4: TTN data for the scaling function Fvn(x). Correc- 
tions to the scaling are negligible and all data collapse. Inset: 
Same data in log-log scale showing the power-law behavior for 
small x with the predicted exponent 1 /4 and the prefactor n. 




FIG. 5: TTN spectrum of the reduced density matrix. In the 
scaling variable of the horizontal axis all data collapse on the 
CFT prediction (7). 



Finally, we consider the full spectrum of pa- If the 



moments of pa behave like Trp^ ~ L 



-e/6(n— 1/n) 
off 



with a 



prefactor roughly independent on n, then the spectrum 
displays the super-universal (i.e. independent on any de- 
tails of the theory) form [25] 



n(A) = / dXP(X) = 7o(2Vbln(A m /A)) 



(7) 



where n(X) is the mean number of eigenvalues larger than 
A, A m the maximum eigenvalue, b = — lnA m , and Io(y) 
a Bessel function. This implies that if n(X) is plotted 
against y = 2-y/61n(A m /A) all data of any system should 
collapse on the same curve. In Fig. 5 we plot n(X) against 
y and all TTN data at different L, £, r (for a total of more 
than 10 5 points) collapse on the curve predicted by CFT. 
Finite size effects are present for small I. Such good 
agreement is due to the fact that (? n F n {x) slightly de- 
pends on n, varying by few per cents in the range [2, oo]. 
This spectrum is fundamental to describe the scaling of 
numerical algorithms [26]. 

To summarize, we reported a full analytic and numer- 
ical analysis of the entanglement of two disjoint intervals 
in the Ising universality class. This represents the first 
numerical check of the CFT predictions (also derived in 
this letter) for quantities that are more complicated than 
the entanglement of the single block. It would be inter- 
esting to understand how these results change in sys- 
tems with boundaries (that already for the single inter- 
val present intriguing features [9, 27]) and in the pres- 
ence of quenched disorder, to understand if the apparent 
"restoration" of conformal invariance for one interval [28] 
is somehow preserved in the case of many. 
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